On "jamitons," self-sustained nonlinear traffic waves 
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Abstract 

"Phantom jams," traffic blockages that arise without apparent cause, 
have long frustrated transportation scientists. Herein, we draw a novel 
homology between phantom jams and a related class of self-sustained 
transonic waves, namely detonations. Through this analogy, we de- 
scribe the jam structure; favorable agreement with reported measure- 
ments from congested highways is observed. Complementary numerical 
simulations offer insights into the jams' development. Our results iden- 
tify conditions likely to result in a dangerous concentration of vehicles 
and thereby lend guidance in traffic control and roadway design. 



Continuum traffic models whereby individual vehicles are modeled col- 
lectively have enjoyed a rich history since their popularization by M. J. 
Lighthill and G. B. Whitham. These models are able to explain dynamics 
both on open roadways and in the vicinity of traffic lights and also elucidate 
"phantom jams," blockages that arise without obstructions [1. Notwith- 
standing this venerated record, the possibility of self-sustained disturbances 
consisting of a shock matched to a transonic flow has received inadequate at- 
tention despite recent experimental evidence that such "jamitons" are ubiq- 
uitous [2J. We demonstrate herein that jamitons are structurally similar to 
detonation waves [3]. This insight, fascinating in its own right, lends helpful 
guidance in roadway design. 




Figure 1: Profiles of the non-dimensional traffic density for 
[p-/pM, «-/«o) = (0.1, 0.9), (0.2, 0.8), ... (0.5, 0.5), where p_ and 
U- denote the traffic density and speed far up- or downstream of the shock. 
Here p x /p = (3p x /(p M - p) and u = u (l - p/pu) with (3 = 10m 2 /s 2 , 
Pm = 0.2m" 1 , uq = 20m/s and r = 5s. Stars denote the sonic point. 



Continuum traffic models typically consist of coupled pair of partial dif- 
ferential equations describing the mass and momentum of the flow |T] . Using 
subscripts to indicate differentiation in time and space, this coupled pair is 
expressed, generically, as 

Pt + {pu) x = , (1) 
u t + uu x + — = -(u-u) , (2) 

P T 

where r is a relaxation timescale and u and p are, respectively, the local 
traffic speed and density. From the discussion of Aw & Rascle [3] , the traffic 
pressure, p, increases with density. The desired speed, u, decreases with 
density and vanishes as p approaches the maximal density, pu- 

In the experiment of Sugiyama et al. [2] , 22 vehicles were placed on closed 
0.23 km long track; the distance between adjacent vehicles was initially con- 
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Figure 2: Multiple jamitons. Orange and blue curves show, respectively, the 
theoretical and numerical solutions. The difference is visible only in the ex- 
panded inset, which also presents a representative distribution of individual 
vehicles along the bottom. 



stant. Instructing drivers to travel at 30km/h, researchers soon discerned 
spatial inhomogeneities in p and u. Most notably, a propagating density 
bulge, approximately six vehicles thick, was observed to travel against the 
direction of traffic flow. Away from this jam, vehicles moved unencumbered. 
Each vehicle decelerated as it entered the jam, was momentarily delayed, 
then accelerated out the other side. 

Motivated by these observations, self-sustained solutions are sought (on 
an infinite road) in terms of the self-similar variable rj = (x — st)/r where 
s (< u) is the traveling wave speed. In terms of rj, ([I]) and ^ may be 
combined to yield 

du (u — s)(u — u) 
dr] (u — s) 2 — c 2 

Here c = (dp /dp) 1 / 2 is analogous to the sound speed from compressible 
flow. The solution of ([3| follows from analyses in disciplines as distinct as 
combustion [3], astrophysics [5], and arterial blood flow jB]. A sonic point 
is an event horizon where, in the frame of the traveling wave, the vehicle 
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Figure 3: Vehicle trajectories for 
(p-/pM, u-/uq) = (0.35, 0.65). 



Figure 4: Measured vehicle trajecto- 
ries due to J. Treiterer. The green line 
shows the experimental measurement 
of [2] (Figure reproduced from |2J with 
permission) . 



speed matches the speed of small disturbances, and is defined by u — s = c. 
The speed and structure of the jamiton is obtained by integration of Q 
under the condition that du/drj remain finite at the sonic point where both 
the numerator and denominator of ^ simultaneously vanish. The jamiton 
is a self-sustained wave because the region between the shock and the sonic 
point is decoupled from the flow downstream of the sonic point in direct 
analogy with self-sustained ( Chapman- Jouguet) detonation waves [3]. 

Solutions are presented in figure [TJ which shows profiles of the traffic 
density. Consistent with Sugiyama et al. [2], jamitons connect up- and 
downstream states of equal density and speed and arise only when, as in 
figure [TJ the background density exceeds a critical threshold. There is a 
striking correspondence between the curves of figure [TJ and the detonation 
wave depicted in figure 2.1 a of Fickett & Davis [3]. 

As with the roll waves observed in gutters on rainy days, numerical solu- 
tion^] of |l| and ([2]) shows the development of a train of jamitons. Figure [2] 
indicates that individual waves from this train are well-described by the 
foregoing theoretical analysis even when, as with the left-most jamiton, the 
waves have not completely saturated. 

Vehicle trajectories are shown in figure [3] Left-to-right moving vehicles 
travel at constant speed U- until they encounter the upstream-propagating 

1 Details of the numerical method are briefly reviewed in Appendix A. Further detail 
will be presented in a forthcoming publication. Additional examples of numerical output 
are given at http://math.mit.edu/projects/traffic. 
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jam. Thereafter, a rapid deceleration, followed by a gradual acceleration, is 
observed. Data from congested highways show similar trends (figure 

Descriptions of jamitons provide helpful design guidance as they iden- 
tify conditions likely to produce dangerous vehicle concentrations. Such 
situations may be avoided by judicious selection of speed limits, carrying 
capacities, etc. 
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A Numerical method 

To investigate the time evolution of density and speed anomalies within a 
stream of vehicles, the governing equations (1) and (2) are solved numeri- 
cally using a mesh-free, Lagrangian particle method [7] calibrated to accom- 
modate rapid changes in p and u in the vicinity of (numerically-smeared) 
shocks. The continuum of vehicles is represented as a set of O(10 3 ) — O(10 4 ) 
discrete particles whose motions are governed by the combination of forces 
given on the right-hand-side of (2). Individual particles neither appear or 
disappear; consequently the mass continuity equation (1) is satisfied au- 
tomatically. The number of particles exceeds by two or more orders of 
magnitude the number of vehicles. Therefore, although there is a natural 
connection between the particle and vehicle density, the numerical method 
is fundamentally "macroscopic" in that the reaction of individual drivers is 
not explicitly modeled. 
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